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Abstract 

We propose a method for deriving networks from one-dimensional binary cellular automata. The derived 
networks are usually directed and have structural properties corresponding to the dynamical behaviors 
of their cellular automata. Network parameters, particularly the efficiency and the degree distribution, 
show that the dependence of efficiency on the grid size is characteristic and can be used to classify cellular 
automata and that derived networks exhibit various degree distributions. In particular, a class IV rule 
of Wolfram's classification produces a network having a scale-free distribution. 
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1. Introduction 

Cellular automata (CA) have been used to study the critical phenomena of complex systems. S. 
Wolfram has systematically investigated the dynamical behavior of one-dimensional automata and iden- 
tified the following four essential types: homogeneous (class I), periodic (class II), chaotic (class III), and 
complex (class IV) [1]. In particular, class IV rules produce complex structures with long transients. 
Alternatively, complex networks exhibiting a scale-free topology can be identified ubiquitously; for exam- 
ple, they occur in social relationships [2], biological and chemical systems (1,131, and the Internet (H, @]. 
Such networks have been studied extensively in the wake of papers by Watts and Strogatz on small- world 
networks [7] and by Barabasi and Albert on scale- free networks [8j. Even though both CA and complex 
networks are used to study complex systems and various phenomena, their relation is not clear. 

In recent years, the dynamics of Boolean networks with complex topology have been studied j9[- 13| 
yielding several notable results such as the dependence of dynamical phases on network topology 10] 
and the robustness of scale- free networks [11]. Our approach is to define networks that correspond to the 
dynamical behaviors of CA. 

In this article, we propose a method for deriving networks from the time evolution of CA config- 
urations. For transformations of the CA rule function, the adjacency matrix of the network exhibits 
characteristic properties. We investigate networks derived from the typical rules of one-dimensional bi- 
nary CA with three and five neighbors. Our studies reveal that the structural properties of the derived 
network reflect the dynamical behaviors of the CA rule and that chaotic or critical rules lead to complex 
network topologies. We use two parameters to characterize the topology of such networks: the efficiency 
of the network (corresponding to the harmonic mean of its shortest path lengths) and its degree distribu- 
tion. The efficiency shows a characteristic scale dependence on the grid size of the cellular automaton and 
may be useful in classifying CA rules. Class III rules correspond to random networks, and a scale-free 
degree distribution has been obtained from the class IV rule. 

The next section describes our notation and some definitions relevant to CA. In section 3, we describe 
our method for deriving networks from CA rules and discuss their properties under transformations of 
the CA rule function. Section 4 reports on the efficiency and degree distribution of derived networks and 
discusses the correspondence between the dynamical behavior of a cellular automaton and the structural 
properties of its derived network. 

2. Notation and definitions relevant to cellular automata 

CA are dynamical systems that consist of a regular grid of cells, each characterized by a finite number 
of states. CA are updated synchronously in discrete time steps according to a local rule (CA rule) that 
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is identical at every cell. In a one-dimensional grid, each cell is connected to its r local neighbors on 
either side, where r is referred to as the radius. Thus, each cell has 2r + 1 neighbors, including itself. 
The state of a cell for the next time step is determined from the current states of the neighboring cells: 
Xi(t + 1) = f(xi- r (t), ...,Xi(t), ...,Xi+ r (i)), where xi{t) denotes the state of cell i at time t, and / is the 
transition rule function. The term configuration refers to an assignment of states to all the cells for a 
given time step; a configuration is denoted by x(t) = (xo(t), xi(t), Xjv-i(t)), where N is the grid size. 
Thus, the time transition of a configuration x(t) with periodic boundary conditions is given by 



x(t + i) = f(x(t)) (i) 

= (f(x N - r (t), . . . , X (t), . . • , X r (t)), 

f(XN-r+l(t), • • • , Xi(t), . . . , Xi+ r (t)), 

. . . , f(x N -i- r (t), . . . , X N -i(t), . . . , X r -l(t))), (2) 

where / represents a mapping on the configuration space {x} N . In this article, we restrict our 
discussion to binary CA, which satisfy X{ G {0, 1} for all i. 

For a given configuration a?, the mirror (left-right reflection) and complement (0-1 exchange) config- 
urations are denoted as x = (xjv-i, xi, xq) and x = (xo, xi, #jv-i)> respectively. By analogy, the 
mirror, complement and mirror-complement of a rule function / are defined as follows: 



f (Xi— r i - - - 1 Xi , . . . , Xi-\- r ) — f (xi-^-rp , . . . , Xi , . . . , Xi — r ) (3) 
f(Xi— r 5 •••■) Xi) 3^-|_ r ) = f(Xi— r , . . . , Xi , . . . , X^-j-r) (4) 



f(Xi— r 5 •••■) Xi) 3^-|_ r ) — f(Xi+ r , X^, X^_ r ) , (5) 

respectively. The CA rules of these transformed functions are equivalent to the original rule / [14] . It is 
trivial that mirror and complement operations are commutative, i.e. / = /. The mappings defined from 
these functions are 

f(x) = (f(x N - r , . . . ,X , . . . ,X r ), f(x N - r +i, . . . ,Xi, . . . ,Xi +r ), 

. . . , /(xAr_i_ r , . . . , arjv-i, • • • , x r _i)) (6) 

= (f(x r , . . . , X , . . . , XAT_ r ), f{x\+ r , . . . , Xi, . . . , XjV-r+l), 

. . . , f(x r -i, . . . , XjV-1, . . . , XjV-l-r)) (7) 

= /(*) (8) 

/ 0*0 = (f(xN-r, • • • , X , . . • , X r ), f(x N - r +i, . . . , Xi, . . . , X i+r ), 

. . . , f(x N -i- r , . . . , Xjv-l, . . • , X r _i)) (9) 
= (f(x N - r , . . . , X , . . . , X r ), f(x N - r +i, . . . , Xi, . . . , X i+r ), 

. . . , f(xN-l-r, • • • , • • • , X r _i)) (10) 

= (ii) 

/(*) = m = }(x). (12) 

A £-fold repetition of these mappings yields, respectively, 



/'(*) = (13) 
f\x) = fU*l (14) 
}\x) = fUx) = r(x). (15) 

Elementary Cellular Automata (ECA) are the simplest nontrivial binary CA; they are defined on 
a one-dimensional grid with minimal neighborhood size (r = 1). The 2 3 = 8 different neighborhood 
configurations result in 2 8 = 256 possible rules, of which 88 are nonequivalent under the transformations 
©-(O ECA rules are generally referred to by their Wolfram code, a standard naming convention 
invented by Wolfram P, [TEJ that gives each rule a number from to 255. For example, rule 30 exhibits 
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(a) rule 62 (b) rule 30 (c) rule 90 (d) rule 54 (e) rule 110 

Figure 1: Examples of network graphs derived from EC A rules with N = 41 and t = 20. The dots on the 
circumference of each graph correspond to the nodes of the derived network, which in turn correspond 
to the cells of the cellular automaton. The time t is selected to give each cell causal relationships with 
all other cells and to avoid repetitions. 



class III behavior, meaning that even simple input patterns lead to chaotic, seemingly random histories. 
Rule 90 is also chaotic with a fractal structure (Sierpinski triangle). Furthermore, rule 110 generates class 
IV behavior, which is neither completely random nor completely repetitive. Localized structures appear 
and interact in various complicated ways. 

Other simple models are 5-neighbor (r = 2) CA, which contain 2 32 rules. Our discussion is restricted 
to the totalistic CA (5TCA), in which the state of each cell at time t depends only on the sum of the states 
of the cells in its neighborhood at the previous time. Of the 2 6 = 64 totalistic rules, 36 are independent. 
To avoid confusion, we add the letter " T" to the Wolfram code of the 5TCA rules, e.g., rule T20. 



3. Derivation of networks from CA rules 



We consider a one-dimensional grid with N cells, where each cell is connected to its r nearest neighbors 
with periodic boundary conditions. Each cell state x G {0, 1} evolves by an identical CA rule function 
fn, where R denotes its Wolfram code. After t time steps, the configuration of cells obtained from an 
initial one cp = x(0) is given by x(t, cp) = If <fi denotes the initial configuration with a changed 

state for cell \, the difference between the configurations after t time steps from (p i and cp can be written 
as 



Aixfa cp) = x(t, cp,) + x (t, cp) (mod 2) = f R (cp z ) + f R (cp) (mod 2) 



(16) 



where Aix(t, ip) characterizes the influence of cell i on other cells in the grid after t time steps. This 
influence represents the flow of information from cell i through the network. We define a matrix as 



^fl (*,¥>) = [A x(t,<p),A 1 x(t,ip),...,A N - 1 x(t,ip)} J 



(17) 
(18) 



where the transpose operation T does not apply to the individual elements Aix(t,(p). We treat 
<p) as the adjacency matrix of a network derived from the CA rule R; an entry = 1 if a 
directed edge from node i to node j exists, and otherwise. Because Aix(t, cp) depends on the initial 
configuration, An(t,cp) gives a different network for each initial configuration. Although this derived 
network is not identical with an analogous network derived by considering the change of two or more 
cells, it captures essential properties of its cellular automaton, as shown below. Some graphs derived 
from EC A and 5TC A rules are presented in FigJT] and Fig]2j respectively. 

From Eqs. ([T3lh ([16]) and ([T5]h the adjacency matrix is transformed by the mirror operation as follows: 



A R (t,cp) 



A /flfoO, Ai/ateO, • • • , A N _J R (cp) 
Ar&&) : 



(19) 
(20) 
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(a) rule T40 (b) rule T10 (c) rule T42 (d) rule T20 (e) rule T52 

Figure 2: Examples of network graphs derived from 5TCA rules with N = 41 and t = 10. 

where the subscript of A in Eq. ([T9)) labels the state-changed elements of tp and the ~ operation acting 
on the adjacency matrix in the right-hand side of Eq. ([20|) represents the mirroring of elements about 
both the horizontal and vertical axes. Similarly, the complement of the adjacency matrix is obtained 
from Eqs.fjUj), (USD and JTBD as follows: 



A R (t,<p) 



A f t R (cp),A 1 ? R (v),...,A N _ 1 f t R (cp) 



Ao/^),Ai^(^),...,Aiv-i^(^) 

= [Ao/^)^!/^),...^^/^)]" (21) 

= A R (t,(p), (22) 

where Eq. (|2T|) is derived from the invariance of the modulo operation in Eq. ([T6|) under 
complementation. In general, the graph of the complement A R (t, cp) is not identical with the "graph 
complement" [16| of A R (t, (p). Finally, the mirror-complement of the adjacency matrix is 



A R (t,<p) = A R (t,$) = A R (t,<p). (23) 

An important property of the adjacency matrix may be derived for self-complementary rule functions. 
Adding the transformations (|3])-(|5]), we define the diminished-radix complement of the rule function f R 

as 



f R (Xi- 



. . . , X{ , . . . , Xi-\- r 



f R (Xi- 



. . . , X{ , . . . , Xi-\- r 



(24) 



If a mapping f R is self- complementary, i.e. f R — f R , the mapping f R obtained from f R satisfies the 
equation 




lor even t. 

Therefore, Aif R (<p) is equal to Aif^tp) for all t, and the adjacency matrix A R (t,(p) defined by f R (<p) 
is identical with A R (t, cp). Although this property of the adjacency matrix makes the rule f R 
indistinguishable from rule i?, this degeneracy is not a drawback of our method but rather a new way of 
detecting similarity among CA rules, as described below. 

Another interesting property of our approach pertains to additive mappings which satisfy the equation 
f R (x + y (mod 2)) = f R (x) + f R (y) (mod 2). Eq.(l6]) leads to the result Aix(t, (p) = f R (0i), where 0^ 
is the null configuration except at cell i. Thus, the adjacency matrix JTHJ) is independent of the initial 
configuration and all nodes of the derived network are equivalent. Hence, the same number of edges 
connects to each node. If an edge from node % to node j exists, an edge also exists in the opposite direction; 
the network is undirected, which means that the adjacency matrix is symmetric. For example, rule 90 is 
additive and, for this reason, corresponds to the geometrical graph shown in Fig. 1(c). Furthermore, if cell 
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(a) ECA pairs 
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Table 1: The diminished-radix complements pairs of self-complement rules of ECA and 5TCA. The rule 
inside the parentheses is mirror equivalent to the one outside the parentheses. 



i and cell j change their states in the initial configuration (p, the difference between the configurations 
<Pij and cp after t time steps satisfies the equation 



A ijX (t, <p) = fniiptj) + f R {<p) (mod 2) = /^(Oy) = f R (0i) + f R (0j) (mod 2) 

= Aix(t, 0) + Ajx(t, 0) (mod 2) (26) 
= (A R (t))i + (A R (t))j (mod 2). (27) 

where (A)i denotes i-th row vector of matrix A. Consequently, the adjacency matrix for additive 
mappings can represent the influence of changes at more than two cells. 



4. Properties of derived networks 

Following the Wolfram's classification, we examine some typical rules of ECA and 5TCA. There is 
no reason to restrict our discussion to even- numbered rules. We use two basic parameters, the efficiency 
and the degree distribution, to investigate the properties and complexities of derived networks. 

A descriptive network parameter is the distribution of degrees. Pr(^ u ) and PR(k out ) indicate the 
probabilities of a node having an in-degree ki n and an out-degree k outl respectively. In particular, scale- 
free networks are the class of networks whose degree distribution is a power-law: P(k) ~ /c -7 , where 7 
is called the scale-free exponent. Small-world networks can be categorized by their average shortest path 
length / =< dij >, where d^ is the length of the shortest path between node i and node j. To sidestep 
any divergence of the dp, we consider their harmonic mean, which we use to define the network efficiency 

Because these parameters are calculated for adjacency matrices obtained from randomly generated 
initial configurations, the transformed matrices defined in Eqs. ([20|h ([22)) and ([23)) give assuredly equivalent 
results with the original adjacency matrix. Hence, the dependence of efficiency on grid size N is very 
useful for the classification of CA rules. Further insight is obtained by using the degree distribution to 
characterize the CA rules; for example, random and scale- free distributions are found in networks derived 
from chaotic and complex rules, respectively. 

Pairs of the diminished-radix complement of self-complementary rules of ECA and 5TC A are listed in 
Table 1. Although the individual rules of a pair have different CA patterns, they have similar statistical 
properties and belong to the same Wolfram class. As indicated by Eq. ([25]h this correspondence is obvious 
from the invariance of the time evolution under the exchange of zero and one (complementation). 

Networks derived from ECA and 5TCA rules 

As shown in Figs JT] and [2j the topology of the derived network reflects properties of the corresponding 
CA rule. Fig|3] and FigHl plots the efficiencies of illustrative ECA and 5TCA rules, respectively. The 
efficiency values can be categorized into several types. All class I CA have zero efficiency, whereas class 
III CA are characterized by a high, TV-independent efficiency (i.e. short average path length), except 
for rule 90, 60 and T42. For class II CA, the efficiency decreases as A/ --1 , with the exception of rule 
184. This TV-dependence of class II efficiency can be explained intuitively as follows. Because nontrivial 
patterns of class II CA are localized, the derived networks are disconnected; hence, the number of edges 
is proportional to N, whereas the total number of node pairs is proportional to N 2 . Therefore, the 
contribution of edges to the efficiency decreases as N~ x . Although it is somewhat difficult to decide 
empirically how to classify rule 41, the efficiency values above N > 1600 provide evidence that it belongs 
to class II. 
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Figure 3: Efficiencies of networks derived from EC A rules, representing the average of ten sampled 
networks. 
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Figure 4: Efficiencies of networks derived from 5TCA rules, representing the average of ten sampled 
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(a) rule 62 with N = 6401 and t = 3200 (b) rule T13 with N = 6401 and t = 1600 

Figure 5: In-degree distributions of networks derived from EC A and 5TCA class II rules. 
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Figure 6: In-degree distributions of networks derived from ECA and 5TCA class III rules. 



The efficiency of rule 110 and T52, which is thought to belong to class IV, varies as TV -0117 and 
iV -0,53 (as estimated by least-squares fitting), respectively; these TV-dependence indicates that they lies 
on the "edge of chaos" [19|. In contrast, rule 54 appears to be indistinguishable from class III rules. 
Furthermore, the efficiency of rule T20 is almost proportional to TV -1 for TV>1600. Thus it is difficult to 
distinguish rule T20 from class II rules using only its efficiency dependence. 

Rules 90, 60 and T42 are exceptional. Because they are additive, each node has the same number 
of edges, and the total number of edges is proportional to N. Hence, their efficiencies have similar 
tendencies to those of class II CA. However, in the case of t = 2 n — 1 (where n is a positive integer) 
and TV = 2t + 1, each node of rule 90 has 2 n = (TV — l)/2 edges, which means that its numbers of total 
edges are proportional to A^ 2 . Consequently, the efficiency of rule 90 oscillate between class II and class 
III regions. Similar oscillations of the efficiencies of rule 60 and T42 are confirmed. Another exception 
is rule 184, which is usually called the "traffic rule"; its derived network has random values of efficiency, 
indicating strong dependence on the initial configuration. This random behavior can be interpreted as 
critical phenomena of the phase transition described in Wolfram [15j . 

The in-degree distributions of rule 62 and 7T3 (class II) and those of rule 30 and 7T0 (class III) are 
illustrated in Fig|5]and Figj6l respectively. Rule 62 is a class II rule but has a large components of class 
III character [20], which may account for its long-tail distribution. Rule 30 is well-known for its random 
behavior, from which a random network with a Poisson distribution of degrees is derived. Similarly, rule 
7T0 has two peaks, each fitting a Poisson distribution; this double-peak property corresponds to the 
double-triangle structure of its CA patterns. 

FigJTl shows that class IV rules exhibit exponential decays and long-tail distributions. Figj8] is an 
unaveraged distribution of ten sampled networks of rule T20, from which we obtain the scale- free exponent 
7 = 2.096, based on linear fitting. Hence, a scale-free network can be derived from a class IV rule. 
However, the N~ x efficiency dependence of rule T20 suggests that the network derived from rule T20 is 
disconnected in many components, which contrasts with other scale-free networks such as the BA model. 
We have confirmed that the network derived from rule T20 with = 6401 and t = 1600 has dozens 
of disconnected components. In other words, the network can be rendered scale- free and connected by 
adding, at most, a few dozen edges. By examining Figs|6{b), W[b) and [T^a), it is apparent that the 
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(a) rule T20 



(b) rule T52 



Figure 7: In-degree distributions of networks derived from 5TCA class IV rules with N = 6401 and 
t = 1600. (a) The distribution of rule T20 plotted on log-log scale exhibits scale-free behavior. 
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two peaks gradually collapse; the left peak becomes the scale-free distribution, whereas the right peak 
becomes a structure at the lower slope. This change of network structure may represent a phase transition 
from a chaotic phase to a periodic phase. 

5. Conclusions and discussion 

In this article, we have proposed a method for deriving networks from binary CA rules. Networks of 
representative ECA and 5TCA rules have been shown and we have discussed their properties in terms 
of two parameters; their efficiency and their degree distribution. The efficiency parameters appear to be 
useful in classifying CA rules. Representative degree distributions of complex networks have also been 
determined, exhibiting Poisson, long-tail and scale-free distributions, our approach has other interesting 
aspects, such as a new way to characterize additive rules and the similarity of CA rules. We conclude 
that our derived network is an effective representation of CA. 

A survey of all ECA and 5TCA rules will be reported in a full-length paper. However, our method 
can be applied or extended to other CA rules. For example, it is straightforward to apply it to binary 
CA with two or more dimensions and many neighbors, such as the "Game of Life". The extension of the 
adjacency matrix to CA with three or more states may yield weighted networks. 

One remaining problem is how one is to understand the appearance of a scale-free degree distribution. 
The fundamental elements leading to emergent scale-free properties are thought to be the growth of the 
network and the preferential attachment of its links. The former corresponds to the time evolution of 
CA and the latter may derive from the CA rule and the randomly chosen initial configuration. However, 
if the following correspondences between dynamical phases of complex systems and types of complex 
network are assumed, 

• Periodic phase <^> disconnected and localized network 

• Chaotic phase <=> random network 

• Critical phase <^> scale-free network, 

it may be possible to describe phase-transition phenomena as evolutions and transformations of network 
structures. If so, the scale-free network may be an intermediate structure of the description. In any 
event, these correspondences and the appearance of the scale- free degree distribution warrant further 
investigation. 
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